library(tidyverse)
## Warning: package 'tidyverse' was built under R version 4.0.5
## Warning: replacing previous import 'lifecycle::last_warnings' by
## 'rlang::last_warnings' when loading 'pillar'
## Warning: replacing previous import 'lifecycle::last_warnings' by
## 'rlang::last_warnings' when loading 'tibble'
## Warning: replacing previous import 'lifecycle::last_warnings' by
## 'rlang::last_warnings' when loading 'hms'
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v ggplot2 3.3.5 v purrr 0.3.4
## v tibble 3.1.3 v dplyr 1.0.7
## v tidyr 1.1.3 v stringr 1.4.0
## v readr 2.0.2 v forcats 0.5.1
## Warning: package 'ggplot2' was built under R version 4.0.5
## Warning: package 'tibble' was built under R version 4.0.5
## Warning: package 'tidyr' was built under R version 4.0.4
## Warning: package 'readr' was built under R version 4.0.5
## Warning: package 'purrr' was built under R version 4.0.4
## Warning: package 'dplyr' was built under R version 4.0.5
## Warning: package 'stringr' was built under R version 4.0.4
## Warning: package 'forcats' was built under R version 4.0.5
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
library(plotly)
## Warning: package 'plotly' was built under R version 4.0.5
##
## Attaching package: 'plotly'
## The following object is masked from 'package:ggplot2':
##
## last_plot
## The following object is masked from 'package:stats':
##
## filter
## The following object is masked from 'package:graphics':
##
## layout
library(igraph)
## Warning: package 'igraph' was built under R version 4.0.5
##
## Attaching package: 'igraph'
## The following object is masked from 'package:plotly':
##
## groups
## The following objects are masked from 'package:dplyr':
##
## as_data_frame, groups, union
## The following objects are masked from 'package:purrr':
##
## compose, simplify
## The following object is masked from 'package:tidyr':
##
## crossing
## The following object is masked from 'package:tibble':
##
## as_data_frame
## The following objects are masked from 'package:stats':
##
## decompose, spectrum
## The following object is masked from 'package:base':
##
## union
library(tools)
library(janitor)
##
## Attaching package: 'janitor'
## The following objects are masked from 'package:stats':
##
## chisq.test, fisher.test
library(htmlwidgets)
## Warning: package 'htmlwidgets' was built under R version 4.0.5
library(ggthemes)
## Warning: package 'ggthemes' was built under R version 4.0.4
figure5 = read_csv("figure_5.csv") %>%
mutate(version = ifelse(version == "BA","Barabási Albert", "Empirical")) %>%
mutate(network = str_replace_all(network, "_"," ") %>% str_to_title())
## New names:
## * `` -> ...1
## Rows: 144 Columns: 6-- Column specification --------------------------------------------------------
## Delimiter: ","
## chr (3): version, network, seeding_method
## dbl (3): ...1, N, critical_q
## i Use `spec()` to retrieve the full column specification for this data.
## i Specify the column types or set `show_col_types = FALSE` to quiet this message.
g = figure5 %>%
ggplot(aes(x = N, y = critical_q, colour = seeding_method, label = network)) +
geom_point() +
geom_line() +
theme_classic() +
scale_colour_tableau() +
facet_wrap(~version) +
labs(colour = "Seed Method", y = "Critical q") +
ggtitle("Cascade Threshold Comparison Across Seeding Methods")
ggplotly(g) %>%
layout(width = 800) %>%
saveWidget(file = "nVScascade.html")
## Warning: Specifying width/height in layout() is now deprecated.
## Please specify in ggplotly() or plot_ly()
g %>%
saveRDS(file = "nVScascade.RDS")
figure11 = read_csv("figure_11.csv")
## New names:
## * `` -> ...1
## Rows: 297 Columns: 4-- Column specification --------------------------------------------------------
## Delimiter: ","
## chr (1): seeding_method
## dbl (3): ...1, q, proportion_adopted
## i Use `spec()` to retrieve the full column specification for this data.
## i Specify the column types or set `show_col_types = FALSE` to quiet this message.
g = figure11 %>%
ggplot(aes(x = q, y = proportion_adopted, colour = seeding_method)) +
geom_point() +
geom_line() +
scale_colour_tableau() +
theme_classic() +
labs(colour = "Seed Method", y = "Proportion of Nodes Adopting Innovation") +
ggtitle("20 Step Threshold Propagation Across q")
ggplotly(g)
g %>%
saveRDS(file = "qVSprop.RDS")
library(igraph)
dd <- read.table("..\\mass_data/less_than_1500_nodes/aves-weaver-social/edges.edgelist")
graph <- graph.data.frame(dd, directed=FALSE)
plot(graph)
